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Considering that there are two kinds of three dimensional flows in a 
labyrinth seal, a jet flow and a core flow, theoretical equations are set up 
concerning the motion of each flow. The pressure distribution within the 
labirinth is calculated, when the rotor shaft makes a small displacement from 
the center line of the casing, keeping parallel with it. 

The theoretical values of cross coupled stiffness obtained by integrating 
the pressure under different labyrinth geometries and operating conditions 
through these formulas are compared with the experimental data presented by 
other researchers. The theoretical and experimental results show a 
satisfactory agreement. 


1. INTRODUCTION 

A self-excited rotor whirl sometimes occurs in high speed, high pressure 
turbo-machinery, and it is argued that one of the causes of the vibration is 
due to the labyrinth seal force. Recently, carefully prepared and precisely 
measured experiments have been carried out, which make possible a fairly exact 
prediction of the labyrinth seal forces. However, the scope of validation of 
the prediction should be within the range of reference of the experiments, 
because of the complexity of the phenomena, unless theoretical formulas 
verified with experimental data are available. In reality, commercial turbo- 
-machines are normally furnished with labyrinth seals having a number of mixing 
chambers and are operated under much higher pressure than in experiments. 

This paper proposes a method of calculating the asymmetrical pressure 
distribution in a labyrinth seal caused by the parallel displacement of the 
shaft, introducing a new mathematical model suggesting that there exist two 
kinds of independent three dimensional flows, interacting with each other, and 
that the circumferential variation of static pressure in the mixing chamber 
comes from the lack of uniformity of the circumferential velocity of the core 
flow, which is partly originated from the spiral flow effect of the jet flow. 


2. NOTATION 

ao^co^dz.da.eo = factors relating to fluid friction 
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c , C m/ s 3 » circumferential velocity of core flow 

f r a t C m 2 3 “ cross sectional area of mixing chamber and that of core flow, 
2 respectively 

r s »Cm D = cross sectional area of jet flow on the meridian plane 
between successive seal strips 
N,C-3 = number of mixing chamber 
n , Cm” 1 3 - convergency of seal clearance 
p , CP a 3 = pressure 
p Q , ( P a 3 = inlet pressure 

q » C kg/m-s 1 = mass flow of leakage per unit time per unit 
circumferential length 

q D . C kg/m-s 1= ditto, under concentric position of rotor 
R , C J/kg-K3 = gas constant 

R 0 , CmD = inner radius of casing 
r,CnO - radius of rotor 
s , Cm3 — pitch of seal strips 
s 0 » Cm3 = wetted perimeter of core flow 
T, C K3 = gas temperature 
u , C m/s 1 = circumferential velocity of jet flow 
u D , Cm/s 3 = tangential velocity of entry swirl 
u-p » Cm/s 1 ,as circumferential velocity of rotor 
a,C-\3 - factor of contraction of, flow 
/? , C — 3 = pressure ratio, exit to inlet 
<S,Cm3 = tip clearance of seal strips 
S Q , Cm3 - ditto, in the middle of labyrinth 
0 , C rad } = helical angle of jet flow 
d Q , Cradl= angle of expansion of jet flow in meridian plane 
a 1 , 1 2 ,a, 2 > 1 3i C-D = coeficients of fluid friction in circumferential 
direction for r„r 2 ,r J2 ,r 8 , respectively 
ft, CP a *s 3 “ viscosity of gas 

2 / 1-ia , C— 3 = coeficient of carry over from (i-1 )-th to i-th chamber 

9 C kg/m 3 3 = density of gas 

r i» r 2 » r i 2* r 3 C P a^ * circumferential component of shear stress due to fluid 

friction between jet jEIow and casing wall* between the core 
flow and the rotor, between the jet flow and the core flow, 
and between the jet flow and the seal strip , respectively 
Subscript , Superscript 

i = quantity at the i-th chamber in the direction cp 
1-1, i = quantity at the boundary between (i-l)-th and i-th chamber 
in the direction <p 

= quantity under concentric position of rotor 

3. FLOW OF GAS IN A LABYLINTH 
3, 1 Jet flow and core flow 

It has been known for many years that the meridian flow in a labyrinth is 
somewhat similar to the flow ilustrated in Fig. 1. The flow of gas which has 
passed through the tip clearance of the seal strip, here called a jet flow, 
expands and increases its width when it goes through the mixing chamber. This 
flow also has a circumferential velocity due to the entry swirl and the 
peripheral velocity of the rotor. Thus, gas particles in this flow move 
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downwards along a helical stream line, illustrated in Fig. 2. On the other hand, 
the gas flow in the mixing chamber, here called a core flow, forms a vortex 
within the surrounding walls on a meridian plane. This core flow also moves in 
a circumferencial direction induced by the tangential speed of the rotor, and 
consequently the particles in this flow move along a helical stream line in the 
mixing chamber, as illustrated in Fig. 2. 

3.2 Assumptions on the flow 

The following assumptions are used to analyze the three dimensional 
flows : 

(a) Fluid which flows in the labyrinth is an ideal gas* 

(b) Since the change in temperature caused by a change of pressure 

is neutralized immediately, the temperature in the labyrinth is constant. 

(c) The static pressure is constant on a meridian plane in a mixing chamber. 

(d) The coefficient of "carry over", the coefficient of fluid friction and 
the wetted perimeter of the channel when the rotor is in a eccentric 
position, are to be the same value as those under the concentric 
position of the rotor. 

(e) Pressure difference between two adjoining points partitioned with 
a seal strip is small. 

(f) The boundary of the jet flow in a chamber on a meridian plane is 
a tangent to the outer edge of the throat. 

(g) Interchanging of fluid mass takes place between the jet flow and 
the core flow. 

(h) The circumferential component of the core flow velocity is 
constant on a meridian plane in a chamber. 

(i) Even in the eccentric position of the rotor, the influence of the 
eccentricity on the jet flow velocity is small enough for 

its variation to be neglected. 


3.3 Fluid friction acting at the boundary of the flow 

The shearing stress r > due to fluid friction, acting at the boundary of a 
turbulent flow which flows in a narrow gap between two parallel planes can be 

expressed by the formula, , 

r=|- A, 0 V 2 , X = 0.133Re 4 

Where Y is the mean velocity of the flow, and Be is Reynold Vs number relating 
to the distance of the planes. 

In order to estimate the force acting between the flows in a labyrinth the 
above relationship is applied to the flows. Here, since, the movement of fluid 
particles in both the jet and the core flow in a circumferential direction is 
mainly discussed, it becomes necessary to know the relationship between the 
shearing stress and the velocity, both in a circumferential direction. This 
relationship can be obtained by considering the relation between the circumfe- 
rential flow velocity and the circumferential component of shearing stress 
acting in the direction of absolute velocity. Thus, in case of a labyrinth with 
rotating seal strips and of the conditions, u<c <u r the shearing stress r„r 2 
and r 3 acting between the jet flow and the casing wall, between the core flow 
and the channel wall, and between the jet flow and the seal strip, respectively 
relating to the circumferential direction, can be expressed as follows: 
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( 1 ) 


, , I - 1 

pA|U 2 , A| =0-133 c 1 + (v/u) 3 s C a 0 qou/(/* v )3 4 

r 2 = 1 ^ 0 X 2 (u r -C J, A 2 = 0 - 133 C 1 + ( b q v)V( u r — c ) z 3 8 C ( U r C ) C©f So/O 3 4 (2) 

r,= \ pM “r - “)* * X * = djX « 

By a similar principle, the shearing stress acting at the boundary between the 
jet flow and the core flow in a circumferential direction leads to the following 
equation. 

T i2 = 4 " pAj 2 (c— u) , X | 2 = d,d 2 A l X 2 (/A^ -^/Xzdz) */l+( 1 ~c 0 ) 2 v 2 / ( u — c ) 2 (3) 


4. FUNDAMENTAL EQUATIONS OF THE FLOW IN A LABYRINTH 

4. 1 Equations at the eccentric position of the rotor 

By the small displacement x, y of the rotating shaft , as shown in Fig. 3, 
a small variation in static pressure will be added to the pressure distribution 
of the concentric labyrinth- When the seal strips are mounted on the rotor and 
the circumferential velocity of the jet flow is smaller than the peripheral 
velocity of the rotor, the fundamental equations dominating the flows are 
derived in the following way in order to calculate such a variation in pres- 
sure. 

(1) Seal clearance 

From the geometrical configuration, comes the equation, 

t,i = < 5oC 1 — ns( 2 i — N— 2 ) /% 3 — .x c© &<p — 7 si n<p ( 4 } 

(2) Mass flow of gas leaking through the seal clearance 

Instead of neglecting the approaching speed of gas to the throat, the effects 
of "carry over" of the labyrinth seal is introduced. Thermodynamic relations 
provide : 

— (pi-, — p|) /(RT) (s) 

(3) The equation of continuity 

The stream lines of the jet flow in a mixing chamber are approximately repre- 
sented by straight lines having gradient of d% , as shown in Fig. 4(a). A and 
B is the inlet and the outlet of a stream tube having a small circumferential 
length rd<p at the center of the mixing chamber, and and q B represent the 
mass flow rate at A and B, respectively. Due to the difference in the mass 

flow rate the mass in the stream tube segment AB is increased by 

rd<pdt in the duration of time dt . A part of the above increase in gas increa- 
ses its own density and the rest of the quantity, qsi moves to the core flow 
just beneath the stream tube, so that 

qsi = <qi-i,i--qi r L+i ) ~ f- T tan qi- t ,i+ qi^+i ) ^ 

On the other hand, as illustrated in Fig. 4 (b) , the core flow segment 
which exists between D and E at the time t is to be shifted circumferentially to 
the position between D T and E 1 at the time t+dt # 

Since the increase of mass in the volume element D f E f against that in the 

volume element DE is equal to the quantity of gas having come from the jet 

flow, the equation for the conservation of mass in a circumferential direction 
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provides the following equation of continuity 


. hL <- 

dt 7 


dtp 


(4) The equation for the conservation of momentum of the core flow in a 
circumferential direction 

As the increase of momentum in a circumferential direction in the volume 
element D f E f against that in the volume element DE, shown in Fig. 4 -(b)., 
consists of the increase in the circumferential momentum of the gas coming from 
the jet flow, of that due to the pressure gradient and the increase of the 
cross-sectional area of the core flow, and of that due to the fluid friction, 
the following equation can be obtained. 


C-II d( pfm)i j dc± j j?c-u)cl 3(pfa)i + (2c— u)j 
pfm t d t Tl Ji chp r 


dci 

dip 


RT 

pir 



(c-u ) 2 — ? (ur* 



= 0 
i 


(7) 


where K? = li*S/4fm , 


4.2 Relationships when the shaft is at the center line of the casing 


(1) Coefficient of "carry over" in a straight-through labyrinth 

The following equation, which is a slightly modified version of Kamo tori f s 

equation (ref . 1 ) is used. 







<5i— 2 , 1 — i • a+-.s tanflo 


( 8 ) 


(2) Equation of conservation of momentum for the jet flow 

As already mentioned, it is assumed in this case that .so th at the positive 
direction of T[ Z z t Z lz X* should be determined, ai<ci<u r 

Z, :in such a direction that the jet flow is decelerated 

Zz • " 11 11 the core flow is accelerated 

27z : " fl the jet flow is accelerated by the core flow 

25 : " 11 n the jet flow is accelerated. 

In the stream tube segment AB, illustrated in Fig. 5, the increase in the 
momentum of the gas in a circumferential direction in the volume element A\B f 
against AB in a small duration of time dt should be equal to the sum of the 
momentum added to this stream tube segment by the (i-l)-th and i-th and core 
flow in a circumferential direction. Thus, the following equation is obtained. 

(ut—uf— |)q 0 <it =* T| , £— v t 4* ^$ t i — 1 tando^sdt ^i,i+ f*,.! tand 0 ^ a dt 19) 


(3) Fundamental equations when the rotor is in a concentric position 
Substituting x~y=0 and B/dt = d/d<p- 0 lot the fundamental equations, give 

2i-i,l = <y o [l“(2i-N-2)] M 

q|= ( a-3 1-1,1* Ofl 

A A A A A . 

<i( ci — ui) — r'i(u r — ci) —0 02 
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5. SOLUTION OF THE EQUATIONS 

5.1 Mass flow rate, pressure, and circumferential velocity 
of flow at the concentric position of the rotor 


These parameters can be calculated in the following way* 

(a) Values of p t and q Q 

From eqs. (8), (10) and (11), qo can be written as follows: 


^2 

2 _ pQ PW+I 
q °" ~ 
i?i Bl 

RT 


Bi = 


CaSoV 


|l-y(2i-N-2)J 2 -^- jl-y(2i-N-4)J 2 A 1 -,( 2 -A i - 1 ) 


03 


Using the pressure at the inlet p = p , and at the outlet 

A . r o 7 

p N+l = PP 0 , values of qo and p~p N can t> e obtained with the 

necessary accuracy from the above equation through iterative 
calculation. 

(b) Values of ui and c\ 

From eq. (12) 

^ ^ C ~ 

ct-A^iUr + A^iUi 

and, from eq.(14)and eq. (9) substituted by eqs.(l)-(3) 


,1 — Ei)l^ 


8 RT 




a* 


( /\ /\ v 2 ^ /v /\ xv ^ 

, spi-, Ei-,-u r \^ su r (piEi+pi-.i -Ei-i)] 

" qo+ liS — V ul - T° 

are derived, where 


Ei = X ! 2 i i{A2 T i } 2 + <$3 *!,i tan 0 Q 
K /s/s o 


, A 2 ,i=l^Ai,i 


A,,i = 


1 + 




+K/s/s 0 


L l u i c i ' 


0$ 


For the given value of parameters at the (i-l)-th chamber, ui 
can be obtained by solving the quadratic equation (15) for ui . 
However, the values and \ l21 which appear in the coefficient 
of that equation, are unfortunately a complicated function of ui , 
so that ui can be obtained from that equation by applying an 
iterative calculation for each chamber. 


5.2 Pressure distribution when the rotor is located 
in the eccentric position 

The variation of p, c and q due to a small amount of displacement of the 
rotating shaft, must be small compared with those values at the concentric 
position of the rotor. Further, each increment of p, c and q must be a 
periodic function of <p , so that the increment can be expanded by Fourier 
series. Taking the terms of lowest order of the series., the solution of the 
fundamental equations (5) (6)^(7) can be simply put as follows. 
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Pi-pi-r pel ai cos9?+bi si n<p) 

A , 

Ct 555 ct+ u r {di cos^7+ ei s in^? } 

qi-i t i = qo +q 0 { k i-i,i cos^+ U~i,i sin<p} 

where unknown parameters a i,bi...are the linear function of x,y ,and 

|*i cos^-f bi sin^P j <Cpi/p 0 

| di 008^+ ei &in(p |<Cci/u r ^ 

|ki-i,ico8^-f- Ii-,,1 sin^|<d % 

Since the equation obtained by substituting eq. (1 6) into eq. (5) must hold true 
for any value of <p , the following equations can be obtained. 

qo k i-l,i = Po D i,l^ Pi-l a i-l "“Pi a i^ Pi— Pi^ D 2,l x 

q* 1 Pi-i b i-., — Pi bi) — (pi-| — Pi)P 2 ,i y ® 

As similar equations can be also obtained from eqs.(6) and (7) both substituted 
by^f-Ojby eliminating k* 1 ,d, e from these six equations obtained, 

the following equations result eventually. 


Ai,i-i ai-i + Bi,i_i bt-! +At,iai + Bi f ibi+Ai,i + iat+i +Bi, 1+ ,bt-i.| ^Xi-x+Yry 

+Ai ?1 ^ I b 1 -., — Bi^ai -f Ai^ibt —81,14. 1 ai+ i +Ai t i-i. l b 1+ | — — Yi -:x+X t -y 

whe re 


Ai t i-i = -( 1/ q 0 )RXr Po D tjl Pt-i 

At, i = 2p 0 r mi f B<i C Ui - ui Ui+RT 3 (u P - ui)/Mi+ ( l/q 0 ) RTr Po ( 6,, 1+1 +Di,i) pi 

A A 

Ai, i+! ~ ( 1 / q 0 ) RT r Po Oi , 14*1 Pi4*l 

Bi,i-i-=JRT6podiPi,i pi— \S C-2q 0 ) 

Bi,i= — p 0 ?m,i C («i — ill) Ci+RT>(2ci — at) /Mi + RTs p 0 JiC Di,i4i — D t ,i3pi/(2q 0 ) 

A A A . 

Bi f i+i — -RTs po JiD,, i+I pi+i/(2q 0 ; 

Xi s= 2(l— a) a rmi pi Ci(ci — ui)(u r — ui)/Mi 

— ( 1/qo) RTr CB 2 ,i pl«i — {D 2 ,i + D 2l i 4 i } pi + D 2 i i-m p| +i 7 
Yi~ ( 1 —a) spiCiC 1 — (ci— ui) (2ci — ui) /Mi 3 

4*RTs Ji C D 2 ,i Pi— i +($2,14*1 “$2,i> Pi “ $2,14-1 Pi+i .3 / (2qo ) 

Mi - ( 2rm|(u r -ui) J + (2ct — ui) 

D, t i=a 2 5o^i- l , i {l-^ (2i-N- 2)}V(RT) 

B*,i = a* d„ SjL|,* { 1 - ^ (2 i — N - 2 )}/(RT) 

fm,l“ s (^Ro — r — ^ • tznd 0 — aS 0 { 1 — ( 2 i — N — 2) J^j 
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^ ^ A /N 

Ji=ui/vi m i = )C , A2 f i+c'A lt i 

Equations (19) can be written for i=1 upto i=N, so that the total number of the 
equation is 2N. Since Cai-,Di=i = Cb t -0 laM = s Ca i+ ,^ 1 * I3 = Cb i+ ,3i=N = 0 
is evident, the unknown quantities are ai-a H and . And the total 

number is also 2N. The solution, at and bi , therefore, can be obtained by 
solving the simultaneous linear equation (19) considering the following rela- 
tionship. 

at— , fcl^&x+aiy 

The static pressure distribution in a labyrinth can be calculated from the 
eq. (16). 


5.3 Cross coupled stiffness 


The pressure in the labyrinth, expressed by eq.(16), provides the forced 
and f 7 9 acting in the negative direction of x and y axis respectively so that 

n r 2 jt d 

Fx= 2 3T f pi cos<pd<p=7t srpojai 

1=1 fcfl 

Fy s-r J** Pi sin^d^— TTsrp^bi 

The forward tangential force divided by the displacement of the rotating shaft, 
so called the coefficient of the cross coupled spring constant K xy , can be 
represented as follows. 

H 

K X y=*- 3 Tsrpo .2^1 


6. THE FUNDAMENTAL EQUATIONS FOR THE LABYRINTH HAVING 
DIFFERENT GEOMETRIES OR OPERATING CONDITIONS 

The preceding theory refers to the case when the seal strips are 
mounted on the rotor and when the peripheral velocity of the rotor is larger 
than the circumferencial velocity of the jet flow. If, on the same geometry of 
the labyrinth, the peripheral velocity of the rotor is smaller than the 
circumferential velocity of the jet flow, the positive direction of r, and r !2 
should be determined to be in such a direction that the jet flow is decelera- 
ted. Similarly, the positive direction of r 2 should be the direction in which 
the core flow is decelerated. The equations can be obtained, in this case, by 
the same procedure as explained in the previous section. Among these equations 
thus obtained, the following equations are the essentially different as 
compared with those of the preceding section. 

= 0.133 [l+(b 0 v)V(c-u r ) 2 j 8 [(c-« r ) c 0 f a p/( 8 0 /f )]‘ 7 ( 2 a) 


- A A 

1 s Pi(*i ,i + Ei) 

]u*lfa ’PiEiUrV , p 

l 8RT 

J " t + t q ° 4RT J “ 1+ L 


8 Pi-i ( Ei— 


8RT 


/V 


)u 


2 
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_/ -u r \ 

V q o 4RT ) 

fc— u'l ^(pfn)! , dc± ± f (o- u)c') dCpC^i . (2c— u) t 3c t 

+ M_ ^ -[t'(u-< ! ) ! -t'(i!-« r ) ! ) l = 0 <Ta) 

For the straight-through labyrinth, in which the seal strips are mounted inside 
the casing, the equations can also be obtained by modifying the original 
equations with the same principle. 


1-1 8RT J 


(14a) 


7. RESULTS OF NUMERICAL CALCULATION AND THE COMPARISON 

WITH THE EXPERIMENTAL DATA PRESENTED BY OTHER RESEARCHERS 

Factors, a o ,b 0 ,c o ,d h d 2 ,d 3l e o which are included in eqs. (1)-(3), as 

well as a and 0 O included in eq. (8) , are to be decided by the observation of 
the flow in a labgrinth. However, as some of this data was not available, these 
factors have been tentatively estimabed as follows throughont the calculation 
based on the procedure here explained. 

ao — 1-5 b o = 0-5 c 0 = l.Q d , = 1.0 d 2 — 1-0 d 2 — 1-0 eo^ 0-5 
Wachter and Benckert (ref . 2) have done an experiment of the static characte- 
ristics of a labyrinth having different geometries under different operating 
conditions. Results of numerical calculation in accordance with the procedure 
explained in this paper have been compared with the results of the above 
experiments. The results of their experiments shown in Fig. 6 and in Fig. 7 (a) 
are obtained using a land-and-groove labyrinth, whereas numerical calculation 
is on a labyrinth which is deemed to be equivalent to this kind of labyrinth 
and assumed complete mixing of gas in the mixing chamber on making calculation. 
As to the effect of the peripheral velocity of the rotor on the lateral force 
excitation constants, as well as the effect of entry swirl on the constants, 
calculations not only agree with experiments qualitatively, but also fairly 
good quantitive agreement is seen between the two. The dimensionless lateral 
force excitation constants k£ taken on the ordinate of these diagrams, and 
the relative admission energy £* taken on the abcissa in Fig. 7 are the 
following values respectively. 

— * __ dpKxy * u ° 

~ r sNp 0 ( 1 -/ 3 ) * ° ~ 2 ( 1 -/? )KT+C 2 RoqoRT/P 0 (Ro- r 2 ) 3 2 

Wachter and Benckert have done a further experiment using a straight- 
through labyrinth with seal strips morticed inside the casing, and examined the 
effect of entry swirl on the lateral force excitation constants. Comparison of 
experiments and calculations have also been made in this case, and shown in 
Fig. 7(b). A qualitative agreement of the two values is obtained. However, from 
a view-point of quantity, the calculated values are much smaller in this case. 

Also, Jenny (ref .3) has presented his theoretical formulas, which are 
composed partly by adopting the empirical relation of parameters. Results of 
calculation on one example of a labyrinth with axial admission using his 
formula and those by the procedure in this paper are compared and shown in 
Fig. 8. These two calculated values agree very well. 
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It is said that, when swirl is applied to a rotating labyrinth in an 
eccentric position, an additional exciting force is generated in it whose 
intensity depends on the intensity of the swirl- Calculation has been made to 
see the influence of the swirl on the rotating labyrinth. The results agree 
with the tendency, as shown in Fig. 9. 

Kurohashi et a/Xref.4) have done experiment to examine the influence of 
the divergency of the seal clearance on the exciting force of a straight 
through labyrinth. The comparison of calculated values with their experimental 
data is illustrated in Fig. 10. With regard to the statical characteristics it 
is seen in both values that the diverging clearance provides a stabilizing 
tendency against the foreward whirl of the shaft, and the calculated value 
agrees qualitatively with the experimental value. 


9. CONCLUSION 

On the assumption that there exist two kinds of three dimensional flow in a 
labyrinth, a jet and a core flow, an analysis of the behavior of the flow was 
made, and a method to calculate the statical behavior of a labyrinth seal was 
shown. The results of calculation on some examples of labyrinths showed a 
fairly good agreement with results of experiments, and it can be said that this 
method is useful to predict the destabilizing force, which originates from a 
labyrinth seal, and acts on the rotating shaft. 
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(a) Jet flow (b) Core flow 
Fig. 4 Flow on the development surface 



Fig. 5 Stream line of jet flow 
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Fig. 6 Lateral force spring coefficient 
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(b) Labyrinth with seal strips inside the casing 

Fig. 7 Lateral force spring coefficient vs relative 
admission-energy of the flow (4/*» 0) 











Fig. 9 Lateral force spring coefficient of a rotating labyrinth 
affected by the imposing of entry swirl 
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(a) Parameters of experiment 



(b) Comparison of experiments and calculations 


Fig. 10 Effect of divergency of seal clearance 

on the pressure distribution in a labyrinth 
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